Prospects for identifying the sources of the Galactic cosmic rays with IceCube 
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We quantitatively address whether IceCube, a kilometer-scale neutrino detector under construc- 
tion at the South Pole, can observe neutrinos pointing back at the accelerators of the Galactic cosmic 
rays. The photon flux from candidate sources identified by the Milagro detector in a survey of the 
TeV sky is consistent with the flux expected from a typical cosmic-ray generating supernova remnant 
interacting with the interstellar medium. We show here that IceCube can provide incontrovertible 
evidence of cosmic-ray acceleration in these sources by detecting neutrinos. We find that the signal 
is optimally identified by specializing to events with energies above 30 TeV where the atmospheric 
neutrino background is low. We conclude that evidence for a correlation between the Milagro and 
IceCube sky maps should be conclusive after several years. 



GAMMA-RAY OBSERVATIONS 

It is believed that Galactic accelerators are powered by 
the conversion of 10 50 erg of energy into particle acceler- 
ation by diffusive shocks associated with young (1,000- 
10,000 year old) supernova remnants expanding into the 
interstellar medium *l|. The cosmic rays will interact 
with atoms in the interstellar medium to produce pions 
that decay into photons and neutrinos. Dense molecu- 
lar clouds, often found in star- forming regions where the 
supernovae explode, are particularly efficient at convert- 
ing protons into pions that decay into "pionic" gamma 
rays and neutrinos. These provide us with indirect but 
additional evidence for cosmic-ray acceleration and, un- 
like the remnants seen alone, there is no electromagnetic 
contribution to the TeV radiation that is difficult to dif- 
ferentiate from the pionic gamma rays. The existence of 
the "knee" tells us that there must exist Galactic cosmic- 
ray sources producing protons with energies of several 
PeV. These "Pevatrons" will produce pionic gamma rays 
whose spectrum extends to several hundred TeV with- 
out cut-off in interactions with the interstellar medium, 
in particular with dense molecular clouds. By straight- 
forward energetics arguments such sources must emerge 
in global sky surveys with the sensitivity of the Mila- 
gro experiment Q. We will argue that one Pevatron, 
MGRO J1908+06, has likely been identified among five 
candidates in the current sky map. 

Supernovae associated with molecular clouds are a 
common feature of associations of thousands of OB stars 
that exist throughout the Galactic plane. Some of the 
first resolved images of sources in TeV gamma rays 
were of the supernova remnants RXJ1713. 7-3946 [H, 0] 
and HESS J1745-290 which illuminate nearby molecular 
clouds to produce a signal of TeV gamma rays Al- 
though not visible to H.E.S.S., possible evidence has been 
accumulating for the production of cosmic rays in the 
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Cygnus region of the Galactic plane from a variety of ex- 
periments [1, @, 0, S, Most intriguing is a Milagro re- 
port of an excess of events from the Cygnus region at the 
10.9cr level 0. The observed flux within a 3° x 3° win- 
dow is 70% of the Crab at the median detected energy of 
12 TeV and is centered on a source previously sighted by 
HEGRA. Such a flux largely exceeds the one reported by 
the HEGRA Collaboration, implying that there could be 
a population of unresolved TeV 7-ray sources within the 
Cygnus OB2 association. In fact, they report a hotspot, 
christened MGRO J2019+37 '2\. A fit to a circular two- 
dimensional Gaussian yields a width of a — (0.32±0.12)°, 
which for a distance of 1.7 kpc suggests a source radius 
of about 9 pc. As the brightest hotspot in the Milagro 
map of the Cygnus region, it represents a flux of 0.5 Crab 
above 12.5 TeV. 

To date, the Milagro collaboration has identified eight 
Galactic sources of high-energy gamma rays. On the 
basis of prior observations some of these sources ap- 
pear to correspond to objects unlikely to be signifi- 
cant sources of the Galactic cosmic rays. For exam- 
ple, three Milagro hotspots are at the same locations 
as the Crab nebula, Geminga, and the Boomerang neb- 
ula. As these objects are known to be pulsar-wind neb- 
ulae, and therefore not likely to be significant proton ac- 
celerators, we do not consider them in the context of 
this study. Three of these sources, MGRO J1908+06, 
MGRO J2019+37, and MGRO J2031+41, have post-trial 
significances of > 4.9 a [l(j] (the only other Milagro 
source of such statistical significance being the Crab neb- 
ula). The remaining two hotspots — candidate sources CI 
(MGRO J2043+36) and C2 (MGRO J2032+37)— are lo- 
cated within the Cygnus region of the Galaxy at Galactic 
longitudes of 77° and 76°, respectively. Another poten- 
tial hotspot, MGRO J1852+01, falls currently somewhat 
below the threshold set by the Milagro Collaboration for 
candidate sources. If confirmed it will be the strongest 
source in Milagro's entire sky map with a flux about 2.5 
times higher than MGRO J2019+37 [Tl|. In the analy- 
sis that follows, we will consider the five identified Mi- 
lagro hotspots as our candidate cosmic-ray accelerators 
and evaluate the impact of MGRO J1852+01 in the event 
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that it is later confirmed as a source. 

We focus in particular on MGRO J1908+06. The 
H.E.S.S. observations of this source reveal a spectrum 
consistent with a E~ 2 dependence from 400 GeV to 
40 TeV without evidence for a cut-off [13] ■ In a follow-up 
analysis 13j the Milagro Collaboration showed that its 
own data are consistent with an extension of the H.E.S.S. 
spectrum to at least 90 TeV (Fig. [1]) . This is suggestive 
of pionic gamma rays from a Pevatron whose cosmic-ray 
beam extends to the 'knee' in the cosmic-ray spectrum at 
PeV energies. Another source with a measured spectrum 
consistent with E~ 2 is MGRO J2031+41 Q- The lower 



flux measured by MAGIC can be attributed to the prob- 
lem of background estimation for Cherenkov telescopes 
in a high density environment like the Cygnus region. 

Not all the sources have known lower-energy counter- 
parts, however. Although the H.E.S.S. telescope array 
discovered a GeV-TeV counterpart to MGRO J1908+06 
and MAGIC to MGRO J2031+41, the VERITAS tele- 
scopes failed to detect an excess at the location of 
MGROJ2019+37 [H. A possible reason for this dis- 
tinction is that this source, located in the Cygnus region 
of the Galaxy, may not be the accelerator but a nearby 
molecular cloud illuminated by a Pevatron beam. While 
the pionic gamma ray spectrum extends to hundreds of 
TeV, it is expected to be suppressed in the TeV search 
window of VERITAS fl6| . Indeed, there could be many 
potential accelerators in the Cygnus region, one of the 
principal star-forming areas of the Galaxy. 

In conclusion, evidence tracing the production of these 
or any other sources of TeV gamma rays to pions pro- 
duced by cosmic-ray accelerators has been elusive. It is 
one of the main missions of neutrino telescopes to pro- 
duce incontrovertible evidence for cosmic-ray production 
by detecting neutrinos associated with the sources. Par- 
ticle physics is sufficient to compute the neutrino fluxes 
associated with the sources discussed. We evaluate in 
detail the sensitivity of IceCube, the first kilometer-scale 
neutrino observatory now half complete, to the Milagro 
sources assuming that they represent the imprint of the 
Galactic cosmic-ray accelerators on the TeV sky. Here, 
we include for the first time in these kind of calculations 
the effect of a finite energy resolution of the detector 
and a zenith-angle dependent angular resolution. While 
the number of events with energies of tens of TeV is rel- 
atively low, we establish that this is optimally the en- 
ergy region where the atmospheric neutrino background 
is suppressed and an excess from these sources can be sta- 
tistically established. While observing individual sources 
may in some cases be challenging, we conclude that evi- 
dence for a correlation between the Milagro and IceCube 
sky maps should be conclusive after several years. 

It is important to emphasize that the photon flux from 
the Milagro sources is consistent with the flux expected 
from a typical cosmic ray generating supernova remnant 
interacting with the interstellar medium (see for instance 
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FIG. 1: The 7-ray and neutrino fluxes from MGRO J1908+06. 
The hollow/shaded regions surrounding the fluxes represent 
the range in the spectra due to statistical and systematic un- 
certainties. Also shown is the flux of atmospheric neutrinos 
at the same zenith angle as the source (dashed line), taking 
into account the source size and angular resolution. 
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FIG. 2: Calculated neutrino fluxes from five Milagro hotspots, 
assuming an E~ 2 flux and gamma-ray cut-off at 300 TeV. 

[H). In other words, the TeV flux is consistent with 
the energetics that are required to power the cosmic- 
ray flux in the Galaxy. Alternative candidates such as 
micro-quasars have been suggested for the sources of the 
Galactic cosmic rays. If that were the case, cosmic-ray 
energetics would require that they leave their imprint on 
the Milagro sky map, but none have so far been observed. 



NEUTRINOS FROM GAMMA-RAY SOURCES 

Determining the flux of neutrinos from measurements 
of a pionic gamma-ray spectrum is straightforward, 
as both are the decay products of pions produced in 
proton- proton collisions. Here we calculate the neu- 
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trino spectra using the method of rjj. It is illus 
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trated in Fig. [T] comparing the gamma-ray spectrum 
from H.E.S.S./MGRO J1908+06 to the calculated neu- 
trino flux at Earth. As the Milagro data extend to 
~ 100 TeV without seeing a cut-off, we take the gamma- 
ray cut-off at 300 TeV, corresponding to a proton cut-off 
at energies of the order of the 'knee'. The calculated neu- 
trino spectra from the five Milagro hotspots considered 
here are shown in Fig. [21 assuming an E~ 2 spectrum nor- 
malized to the Milagro measurement and also assuming 
a 300 TeV gamma-ray cut-off. 

Earlier work on neutrino event rates from Milagro 
sources modeled the proton spectrum in supernova 
remnants and investigated its effect on gamma-ray and 
neutrino fluxes produced inside the accelerators. Given 
the evidence discussed in the previous section, we as- 
sume in this work that the observed gamma rays are not 
produced directly in the sources but in nearby molecu- 
lar clouds resulting in harder spectra which extend up to 
several 100 TeV. 

Neutrino telescopes detect the Cherenkov radiation 
from secondary particles produced in the interactions 
of high-energy neutrinos in highly transparent and well 
shielded deep water or ice with an array of photomulti- 
pliers. They take advantage of the large cross section of 
high-energy neutrinos and the long range of the muons 
produced. The IceCube telescope [19| is under construc- 
tion and will start taking data with a partial array of 
2400 ten inch photomultiplicrs positioned between 1500 
and 2500 m and deployed as beads on 40 strings below 
the geographic South Pole. With the completion of the 
detector by 2010-2011 the instrumented volume will be 
doubled from 0.5 to 1 km 3 . 

The event rate in a detector above a threshold energy 
^thresh from a neutrino flux dN u /dE is given by 



N ( ,. 



= T 



A eff (E u )^-(E„)dE„, 



where the energy-dependent muon-neutrino effective area 
A e ff(E u ) is taken from [2o| . The angular resolution is 
simulated as a function of the zenith angle according 
to [ill and lies between 0.7° and 0.8° for the Milagro 
sources. The energy resolution is assumed to be ±0.3 in 
\og(E„) which seems achievable given the superior per- 
formance of IceCube compared to AMANDA (~ 0.6 in 
log(_E„) taking 21] and accounting for the kinematics 
at the neutrino-muon vertex and the energy losses of the 
muon on its way to the detector). The flux of atmospheric 
neutrinos from the interactions of cosmic-ray protons n 
the Earth's atmosphere, an irreducible background, is 
tabulated in and gives a good parameterization of 
the AMANDA measurements. Also, we assume no sig- 
nificant contribution from the decay of charmed particles. 
We take the size of the search bin to be the radius that 
gives ~ 70% of the measured gamma-ray flux assuming 



MGRO J2019+37 1.4 
MGRO J1908+06 1.1 
MGRO J2031+41 1.6 



MGRO J2043+36 1.5 
MGRO J2032+37 1.5 
(MGRO J1852+01 1.3) 



TABLE I: Angular radius of the IceCube search bin for each 
Milagro source. 
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FIG. 3: Significance of excess above background as a func- 
tion of threshold energy from H.E.S.S./MGRO J1908+06 af- 
ter 10 years. The shaded area represents the uncertainty in 
the H.E.S.S. 7-ray measurements. The Milagro data points 
suggest the lower limit (dashed line). 



Gaussianity of the source emission and the angular reso- 
lution. Table [J lists the search bin radii. Given a mean 
number of background events and a total number of ob- 
served events (obtained using Poisson statistics for the 
sum of signal and background taking into account 30% 
signal reduction), we calculate the probability (p value) 
that the observed number of events is due to random fluc- 
tuations in the background. We define the significance as 
the p value for which 50% of experiments yield an equal 
or lower p value. 

Figure [3] shows the significance as a function of thresh- 
old energy for MGRO J 1908+06 after 10 years' data tak- 
ing. Because of the Milagro data points lying in the up- 
per error range of the fitted gamma spectrum (Fig. [!}, 
the Milagro measurements favor the higher-significance 
range. As the significance of the excess will likely be 
low even after 10 years, it may be necessary to use a 
stacked search that will look for correlations between all 
five Milagro sources of interest and the IceCube sky map 
simultaneously. 

Figure [4^l shows the mean number of signal events in 
IceCube in 10 years from the five Milagro sources (ex- 
cluding MGRO 1852+01) as a function of energy thresh- 
old together with the mean total number of events. The 
significance of the correlation of this catalog with the Ice- 
Cube sky map after 10 years' time is given by Fig. If 
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FIG. 5: Significance of observed excess from the Milagro 
hotspots including the potential hotspot MGRO J1852+01 af- 
ter 5 years of observation time. 



(50%) boosts the significance in Fig.[3J) to 4er (5a). Even 
reducing the energy resolution to ±0.5 in log(E v ) still 
results in a significance of better than 3cr. Furthermore, 
the use of methods such as unbinned searches beyond 
the simple binned method considered here will increase 
IceCube's sensitivity 



CONCLUSIONS 



FIG. 4: (a) Mean number of neutrinos from the Milagro 
hotspots as a function of energy threshold (dashed line) com- 
pared to the background (dotted line) and total mean (solid 
line) number of events from the search bins in 10 years, (b) 
Corresponding significance of observed excess from the Mila- 
gro hotspots. 



the potential hotspot MGROJ1852+01 turns out to be 
real the same significance would be reached after only 5 
years' observation time (Fig. 0. The figures make clear 
that the best prospect for detecting these sources is to 
focus on events above several tens of TeV, where the 
atmospheric background is very low but there are still 
sufficient signal events left. Then, a detection of these 
sources after several years is possible. 

The results we obtain are conservative in several ways. 
The quoted angular resolution is based on simulations as- 
suming AMANDA technology and reconstruction meth- 
ods. A not unrealistic increase in the resolution from 
0.7°-0.8° to 0.5° improves the significance in Fig. UJd to 
4(7. Also, the assumed photon fluxes at the sources might 
be higher due to absorption in the Galactic photon back- 
ground (especially in the Cygnus region), or, in the case 
of MGRO J1908+06, due to currently ambiguous mea- 
surements (see Fig. [T]). An overall flux increase by 20% 



Apart from the Crab nebula, Milagro has clearly iden- 
tified three sources of high-energy gamma-ray emission 
in their skymap. Two of these sources (MGRO J2019+37 
and MGRO J2031+41) are located in the Cygnus region 
and one (MGRO J1908+06) closer to the Galactic cen- 
ter. Furthermore, the Cygnus region contains two candi- 
date hotspots (MGRO J2043+36 and MGRO J2032+37). 
Another potential hotspot (MGRO J1852+01) falls cur- 
rently below the significance threshold set by the Mila- 
gro collaboration. Combining these measurements with 
measurements from Cherenkov telescopes shows that sev- 
eral of these sources have unusual hard spectra consistent 
with E~ 2 , where in one case the spectrum seems to ex- 
tend up to 100 TeV without indication of a cut-off. 

Sources producing such hard spectra extending up to 
100 TeV and more are required to explain the existence 
of the 'knee' in the cosmic-ray spectrum around 3 PeV, 
with young supernova remnants being the best candi- 
dates. However, their observation is difficult as these 
high-energy photons are produced inside the accelerator 
only within the first few hundred years. Indeed, all cur- 
rent flux calculations for the Milagro sources assume a 
gamma-ray production scenario inside the acceleration 
region and in most scenarios predict photon spectra with 
cut-offs below 100 TeV. 

In our paper we adopt the novel idea that a cosmic- 
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ray source can produce a hard gamma-ray spectrum up 
to high energies over a much longer time of several thou- 
sand years if the gamma-rays are produced outside the 
acceleration region in the interaction of the source's Peva- 
tron beam with a nearby molecular cloud [iH ] . Assuming 
an E~ 2 spectrum with a cut-off at 300 TeV (consistent 
with a proton cut-off at the 'knee') we demonstrate that 
IceCube will be able to see these sources after several 
years of observation. For the significance calculations 
we use results from a detailed detector simulation and 
take into account the energy resolution and zenith-angle- 
dependent angular resolution. The former is especially 
important as we demonstrate that the highest sensitivity 
is obtained by specializing to events with energies above 
30 TeV. Similar studies have so far neglected energy mi- 
gration effects. 

With a visibility of 50%, the location of 
MGRO J1908+06 puts it within reach of a future 
km 3 -scale Mediterranean neutrino telescope. However, 
in the southern hemisphere no Pcvatron candidates 
have been discovered so far, perhaps because no all-sky 
instrument like Milagro is currently operational in that 
hemisphere. While a high-resolution pointed telescope 
could resolve what previously appeared to be a diffuse 
source into its individual parts (supernova remnants and 
molecular clouds), it is possible that the high density of 
ambient matter in star-forming regions prevents individ- 
ual sources from dominating the off-source flux to give 
sufficient statistical significance for a pointing telescope. 
In this case a Milagro-like telescope with a broader 
field of view such that background measurements are 
truly "off-source" would be needed. In this context it is 
interesting to note that IceCube may be able to detect 
gamma-rays from the southern sky 2J| and therefore 



be used to search for southern Pevatrons over a broad 
range similar to Milagro's. 
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